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Abstract 

We discuss the transverse momentum Qt distribution of Drell-Yan pair, produced 
in collisions of transversely polarized protons. We calculate the transversely polarized 
Drell-Yan cross section up to 0(a s ) in the dimensional regularization, which gives 
QCD prediction at large Qt- For small Qt, we include all-orders resummation of large 
logarithms due to emission of soft gluons up to next-to-leading logarithmic accuracy. 
At intermediate Qt, the resummation formula is matched with the fixed-order a s 
perturbative results in a systematic way, and we derive the cross section with uniform 
accuracy over the entire range of Qt- 
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Hard processes with polarized beams and/or target enable us to study spin-dependent 
dynamics of QCD. The helicity distribution Aq(x) of quarks inside nucleon has been mea¬ 
sured in polarized DIS experiments, and also AG(x) of gluon has been estimated from the 
scaling violations of them. On the other hand, the transversity distribution Sq(x ), i.e. the 
distribution of transversely polarized quarks inside transversely polarized nucleon, cannot 
be measured in the inclusive DIS due to its chiral-odd nature, 1 ! and remains as the last 
unknown distribution at the leading twist. Transversely polarized Drell-Yan (tDY) process 
is one of the processes where 8q(x) can be measured, and has been undertaken at RHIC-Spin 
experiment. 2 ! 

We develop the QCD prediction of tDY cross section, da / dQ 2 dQ\dyd<f, differential in 
the transverse momentum Qt and rapidity y of the produced lepton pair, as well as in the 
dilepton invariant mass Q and in the azimuthal angle (f> of one of the leptons with respect 
to the incoming nucleon’s spin axis. Although this Qt- and ^/-differential cross section is 
fundamental in view of comparison with experiment, the corresponding formula has been 
unknown so far even in the leading order (LO) in QCD: the lepton-pair production with 
finite Qt via the Drell-Yan mechanism has to be accompanied by the radiation of at least one 
recoiling parton, so the LO term of the cross section is of 0(a s ). The corresponding one-loop 
calculation of the LO term requires the phase space integration separating out the relevant 
transverse degrees of freedom, to extract the cos(2</>) dependence which is characteristic of 
the spin-dependent cross section of tDY. 1 ! j n the dimensional regularization, in particular, 
the relevant phase space integration in D-dimension is rather cumbersome compared with 
unpolarized and longitudinally polarized cases. Furthermore, at small Qt (“edge regions of 
the phase space”), the radiation of soft gluon produces large “recoil logs” that spoil fixed- 
order perturbation theory, and the corresponding logarithmically enhanced contributions 
have to be resummed to all orders in a s . Previous works treated Qt- and/or y-integrated tDY 
cross sections to fixed-order a S) employing the massive gluon scheme 3 ! and the dimensional 
reduction scheme 4 ! ( see a ] so R e fs. 5), 15)). In this Letter, we compute the Qt- and y- 
differential cross section of tDY in the LO (0(a s )) in the dimensional regularization scheme. 
Then the result is consistently combined with the soft gluon resummation of logarithmically 
enhanced contributions for small Qt up to next-to-leading logarithmic (NLL) accuracy, and 
we derive the first complete result of the well-defined tDY differential cross section in the 
MS scheme for all regions of Qt at the “NLL + LO” level. 

We first consider the tDY cross section to 0(a s )\ h/Pi, Si)+/i 2 (-P 2 , -s 2 ) —>► Z(/ci)+/(£; 2 )+X, 
where /q,/i 2 denote nucleons with momentum P \, P 2 and transverse spin Si,s 2 , and Q = 
k\ + k 2 is the 4-momentum of DY pair. Based on QCD factorization, the spin-dependent 
cross section ATda = (do^si, s 2 ) — da(si, — s 2 ))/2 is given as a convolution, ATda = 
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f dxidx2<5H(xi, x 2 ; /i 2 p)ATda(p F ), where ji F is the factorization scale, 5H(x i,x 2 ;/rfi) = 
Si e ii^Qi{ x i] P 2 F )^Qi( x 2 ', fi 2 F ) + <5gi(xi;/ i 2 F )Sqi(x 2 ;/r| 0 ] is the product of transversity distribu¬ 
tions of the two nucleons, and A T dd(y 2 F ) = (da(s 1 , s 2 ; p,%) — dd(s 1 , —s 2 ; y 2 F ))/2 is the corre¬ 
sponding partonic cross section. Note that, at the leading twist level, the gluon does not con¬ 
tribute to the transversely polarized, chiral-odd process, corresponding to hclicity-flip by one 
unit. We compute the one-loop corrections to A F da(iJ, 2 F ), which involve the virtual gluon cor¬ 
rections and the real gluon emission contributions, e.g., q(pi, Si)+g(p 2 , s 2 ) —► l(ki) +T(A; 2 ) + g , 
with pi = XiPi . We regularize the infrared divergence in D — 4 — 2e dimension, and employ 
naive anticommuting 75 which is a usual prescription in the transverse spin channel. 5) In 
the MS scheme, we obtain for the differential cross section 
A F da 


= cos ( 20 )- 


cr 


[X(Q 2 T ,Q 2 ,y) + Y(Q 2 T ,Q 2 ,y)\, 


( 1 ) 


dQ 2 dQ F dyd(j) v ' rj 3N C S Q 2 
with S — (P] + P 2 ) 2 , and the rapidity y = (1/2) ln[(Q° + Q 3 )/(Q° — Q 3 )] where Q 3 is the 
component of the virtual photon momentum Q M along the direction of the colliding beam 
in the nucleon-nucleon CM system. We have decomposed the RHS into two parts: the 
function X contains all terms that are singular as Qt —■► 0, while Y is of 0(a s ) and hnite at 
Qt = 0. Writing X = W°) + as the sum of d(a°) and O(al) contributions, we have 
dC 0) = SH(xi , x 2 ; p 2 F ) d(Q F ), and 
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where /jlr is the renormalization scale, x 3 = \fr e y , x ° = \fr e~ y are the relevant scaling 
variables with r = Q 2 /S, and SPqq\z ) = 22/(1 — z) + + (3/2) 5(1 — z) is the LO transverse 
splitting function. 5 ^ 9) The “+” distribution to regulate the singularity at Qt = 0 is defined 
such that dQ F ([\n n (Q 2 /Qt)\/Qt)+ = ~ (V[ n + 1]) hr n+1 {Q 2 /pt)- For the Y -term of (d, 
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where we used the shorthand notation for the integral that vanishes for Q F = 0 , 
r rl dx 1 
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and defined the variables, yf t q = \[Q^JS + \Jr + Q\j S, xf 2 — x 1,2 \Z^^~Qt7Q^^ x 2 — 
(x\x 2 — x\x 2 )/{;xi — a:] 1 "), following Ref. 6). Note that ^ and x 2 —> x 2 as Qt —> 0, 
and all terms in © are finite when Qt —■► 0 . 

Eq. © with © and © gives the first result for the differential cross section to 0(a s ) in 
the MS scheme. The result is invariant under the replacement y —> —y as it should be. We 
note that, integrating © over Qt, our result coincides with the corresponding Qr-integrated 
cross sections obtained in previous works employing massive gluon scheme 3 * and dimensional 
reduction scheme, 4 -* via the scheme transformation relation 5 * (see also Ref. 15)). For > 0, 
A"* 0 * vanishes. Also, the terms proportional to S(Qt) br ©, including those associated with 
the + distribution ([In n {Q 2 /Qt)\/Qt)+i do not contribute. The cross section © in this case 
is of 0(a s ), and gives the LO QCD prediction of tDY at the large-Qr region: 


A T da LO 

dQ 2 dQ\dyd(f) 


cos ( 20 ) 


a 


3 N C SQ 2 


X^\Ql,Q\y) 


Ql> 0 


+ Y(Q 2 T ,Q 2 ,y) 
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The cross section © becomes very large when Q T -C Q, due to the terms behaving 
~ a s In {Q 2 /Qt)/Qt an d ~ &s/Qt i n ©• Actually, similar large “recoil logs” appear in each 
order of perturbation theory as a™\i\ 2n ~ l {Q 2 /Q\)/Q^, a” In 2n ~ 2 {Q 2 /Qt)/Qt-> an d so on, 
corresponding to LL, NLL, and higher contributions, respectively, 6 *’ 7 ) and the resummation 
of those logarithmically enhanced contributions to all orders is necessary to obtain a well- 
defined, finite prediction for the cross section. We work out this up to the NLL accuracy, 
based on the general formulation 1 * of the Qt resummation. This can be carried out in the 
impact parameter b space, conjugate to the Qt space: the singular term A" of © is modified 
into the corresponding resummed component, which is expressed as the Fourier transform 
back to the Qt space. As a result, the first term in the RHS of © is replaced by 


A r da 


NLL 


dQ 2 dQ\dyd<j) 


= cos (2 0 ) 


a 


3N C S Q 


^e 2 


2 


db—Jo(bQT ) 


xe 


S(b,Q) 


(C qq ® S Ql ) ® 5q z 


b 2 


)Kp) + (*? 


4) 


( 6 ) 


Here Jo(bQr) is a Bessel function, b 0 = 2e 7B with 7 e the Euler constant, and the large 
logarithmic corrections are resummed into the Sudakov factor e' s * b,< ^* with 
S{b,Q) = - f^/ b2 (dn 2 /k 2 ) {A q (a s (n 2 ))ln(Q 2 /k 2 ) + B q (a s (K 2 ))}. The functions A q , B q as 
well as the coefficient functions C qq , C qq , are perturbatively calculable: we get A q (a s ) = 

E“,i (g b,m = E7. (g)" k”>. with 
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at the present accuracy of NLL, and similarly, 

C„(z,a.) = C w (z,a») = 6(1 - z) |l + ^ C F fy-dj + O(aJ) . (8) 

Here we have utilized a relation 8 ) between A g and the DGLAP kernels 5 )’ 9 ) in order to obtain 
Af> of |Q . 10) The other contributions of ©,© have been determined so that the expansion 
of the above formula © in powers of a s (fi 2 R ) reproduces X with @ to 0(a s ). Our results 
© are consistent with the fact that the LL (Hg 1 ' ) ) and NLL (A^q\ B^) contributions are 
universal (process-independent). n ) 

Eq. © with ©, © satisfies the evolution equation controlled by renormalization group/) 
and is formally perfect. As emphasized in Refs. 12), 13) for the unpolarized processes, 
however, several “reorganization” in the formula of this type is necessary for its consistent 
evaluation up to required accuracy, to treat properly the b integration of © that involves 
too short as well as long distance. In the b space, L = ln(Q 2 b 2 /b^) plays role of the relevant 
large logarithm, with b ~ 1/Qt- Expressing the b dependence of the second line in © as the 
corresponding L dependence, and systematically organizing those contributions in the large- 
logarithmic expansion, the result is given by the resummation formula with exponentiation 
of the LL terms a™L n+1 and the NLL terms a™L n in each order of perturbation theory. 12 )’ 13 ) 
However, L becomes large for small b as well as for large b , and thus the resummation 
formula contains the unjustified large logarithmic contributions at large Qt • To reduce 
these unjustified contributions, we use a procedure 12 ) by performing the replacement L —> 
L = In (Q 2 b 2 /b 2 + 1); L = L + 0(l/(Qb) 2 ) for Qb » 1, so L and L are equivalent to organize 
the soft gluon resummation at small Qt , but differ at intermediate and large Qt , avoiding 
the unwanted resummed contributions as L —> 0 for Qb -C 1. 

The b integration in © also receives the contribution from long distance: similarly to 
other all-orders resummation formula, the b integration is suffered from the Landau pole at 
b = (& 0 /Q)e (1 / 2/3oa:s ) Q )) due to the Sudakov factor, and it is necessary to specify a prescription 
to deal with this singularity. 7 )’ 14 ) We follow the method introduced in the joint resumma¬ 
tion: 13 ) decomposing the Bessel function of © into the two Hankel functions, we deform 
the 6-integration contour of these two terms into upper and lower half plane in the complex 
b space, respectively, and obtain the two convergent integrals as |6| —> oo; note, this choice 
of contours is completely equivalent to the original contour, order-by-order in a s , when the 
corresponding formulae are expanded in powers of a s . 

Prescription to define the b integration to avoid the Landau pole is not unique, e.g., 
“6* prescription” to “freeze” effectively the b integration along the real axis is frequently 
used. 7 )’ 10 ) Such ambiguity reflects incompleteness to treat the long-distance contributions 
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in perturbative framework, and should be eventually compensated by the relevant nonper- 
turbative effects. Correspondingly, we make the replacement e s — > e s ^ b ^ F NP (b) in © 
with the “minimal” ansatz for nonperturbative effects 7 !’ 12 !’ 13 ! with a parameter g^p , 

F NP {b) = exp (~g NP b 2 ) . (9) 

We note that the matching of © with the fixed-order result X, which was used to derive 
© and ©, is now violated at intermediate and large Qt due to the replacement L —*• L. 
This can be recovered by subtracting the 0{a s ) contributions that violate the matching. 12 !’ 13 ! 
Combining the result with the LO cross section ©, we obtain our differential cross section 
of tDY at the “NLL + LO” level: 

A T da _ A T dd NLL [ A T dd NLL 
dQ 2 dQ^dyd(j) dQ 2 dQ^dyd(f) dQ 2 dQ^dyd(t) 

where ATdd NLL / (dQ 2 dQpdyd(j)) denotes © with L —■> L and e^’Q) —■> e s ^ b, ^F NP (b), and 
[• • • ]o(a s ) denotes the 0(a s ) terms resulting from the expansion of “• • • ” in powers of at s with 
gNP —> 0. There is no double counting between the resummed and fixed-order components 
in m for all Qt] e.g., at Qt Q , both the second and third terms in m become large 
but cancel with each other. In particular, the integral of m over Qt reproduces that of 
© exactly , because 12 ! L = 0 at b = 0; our Q T resummation does not affect the total cross 
section. Eq. m gives the well-defined tDY differential cross section in the MS scheme over 
the entire range of Qt- 

As an application, we compute (fTUl) as a function of Q T for a /S = 200 GeV, Q — 8 GeV, 
y = 2, and 0 = 0, which corresponds to the tDY Q-r-spectrum for the detection of dimuons 
with the PHENIX detector at RHIC. As the first estimate of this quantity, we use the 
following nonperturbative inputs, for which our knowledge is uncertain: for the transversity 
5q(x,y 2 ) in ©, ©, we saturate the “Soffer bound” as 5q(x,y „) = do) + Aq{x 1 Pq )]/2 at 
low input scale y 0 ~ 0.6 GeV with the NLO density and helicity distributions q(x,y g) and 
Aq(x,y q), and evolve 5q{x 1 yl) to higher /i 2 with the NLO DGLAP kernel, 5 !’ 9 ! see Ref. 15) 
for the detail. As for the nonperturbative parameter g^p of © , we use g^p — 0.5 GeV 2 
suggested by the study of Ref. 16). The solid curve in Fig.l shows (flUj) . multiplied by 2 Qt, 
with g NP = 0.5 GeV 2 . We also show the LO result using © by the dashed curve, and the 
contribution from the Y-term of © by the dotted curve. For calculation of all curves, we 
chose tip = dn = Q- The LO result becomes large and diverges as Qt —> 0, while the “NLL 
+ LO” result is finite and well-behaved over all regions of Qt- The soft gluon resummation 
gives dominant effects around the peak of the solid curve, i.e., at intermediate Qt as well 
as small Qt- In Fig.l, the results of (fTTH) with g^p = 0.3 GeV 2 and 0.8 GeV 2 are also 
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A T da LU 

- dQ 2 dQ‘i,dyd(!) 
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shown by the dot-dashed and two-dot-dashed lines, respectively: the larger g^p gives the 
broader spectrum with the higher peak position, because the larger g^p of @ corresponds 
to the larger “intrinsic transverse momentum” of partons inside nucleon. The impact of the 
nonperturbative effects © becomes milder when the short-distance contributions are more 
dominated, i.e., for the case with larger Q 2 (see also Refs.12), 13)). 

In conclusion the perturbative effects relevant for the dilepton Qp spectrum in the tDY 
in QCD are now under control over the entire range of Qp. Apparently further systematic 
studies in various kinematic region of pp collision, as well as of pp collision, 1 ') with the 
resummation formalism are indispensable to reveal 5q(x) and also the intrinsic transverse 
motion of quarks. 



Q t (GeV) 

Fig. 1. Qp distribution Apda/dQ 2 dQpdydcj) at \/S = 200 GeV, Q = 8 GeV, y = 2 and 0 = 0. 
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